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ABSTRACT 

The order h fluctuations of gauge fields in the vicinity of a blackhole can create a repulsive 
antigravity region extending out beyond the renormalized Schwarzschild horizon. If the 
strength of this repulsive force increases as higher orders in the back-reaction are included, 
the formation of a wormhole-like object could occur. 
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1. Introduction 

The study of wormholes, or handles in the spacetime topology, has been an active research 
area in recent years, with the bulk of the effort going into understanding the behavior of 
topology change in Euclidean space and the putative role these instantons might have for 
solving the cosmological constant problem, spawning baby Universes, and related concerns 
[1]. In physical (i.e., Minkowski-signature) spacetime by contrast, wormholes have attracted 
less attention, probably due to the lack (until recently) of concrete and physically meaningful 
examples. Indeed, the only explicit example for over seventy years has been provided by the 
well-documented Schwarzschild wormhole [2], which is at best a theoretical curiosity. Not 
only does the wormhole throat coincide with the blackhole event horizon, but the "bridge" 
connecting the two asymptotically flat spacetimes snaps into two disconnected pieces so 
rapidly as to prevent any signal (even light) from traversing the throat [3]. However, a 
recent analysis of static isotropic solutions of Einstein's equation revealed these solutions 
contain bona-fide wormholes provided the throat is threaded with stress-energy violating the 
weak energy condition (WEC) [4]. These wormholes are free from the problems plaguing 
the Schwarzschild case. If the requisite stress-energy exists, then one should attempt to pin 
down and expose the physical traits belonging to this class of wormholes and to clarify the 
scenarios leading to their existence. The knowledge derived from this will lead to a better 
understanding of the kind of nontrivial topologies allowed by the gravitational field equations 
and motivates exploring the cosmological consequences these wormholes are bound to have 
in the early Universe. 

A hitch is encountered at the purely classical level however, as no classical stress tensor 
is known to violate the WEC. But the interplay of gravity and quantized matter can lead to 
the type of stress-energy relevant for maintaining a wormhole structure since quantum fields 
in curved backgrounds tend to develop localized negative energy densities [5]. Since the 
production of negative energy densities is a generic consequence of the gravitational interac- 
tion with matter [6] , this lends credence to the picture of the early Universe as a spacetime 
containing naturally occuring wormholes. The task remains to find an explicit wormhole 
solution arising from a physically realized matter source. As a step in this direction, in this 
paper we derive a semiclassical solution of Einstein's equation coming from the influence of 
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quantized gauge fields on the classical background geometry of a Schwarzschild blackhole. 
This means we are to solve the back-reaction problem 



(1.1) 



where < T^j, > represents a renormalized quantum stress tensor for matter in a Schwarzschild 
background (with metric g) , and 5gis the perturbation caused by the source term. The exact 
1-loop renormalized quantum stress tensor has been calculated for U{1) fields [7], and we 
take this as our source term. We also consider the back-reaction of non-abelian gauge fields 
in the limit where the gauge self-coupling is small. In that limit, we can take over the results 
in [7] by noting that the abelian and non-abelian stress tensors are related by 



where C2{G) is the second Casimir invariant (for the group G) evaluated in the adjoint 
representation. If G2{G) ~ 100, (e.g., the string-inspired unification group x has 
C2 = 120) the total effective mass seen by a static observer located within a few times 
the renormalized Schwarzschild radius can be negative, and hence she will feel a repulsive 
gravitational force in that region. As discussed in [4], this repulsion or "defocussing" is a 
necessary condition for the existence of a traversable wormhole. 
2. The Solution of the Back-reaction 

A concise discussion of the back-reaction problem along with details useful for calculat- 
ing on Schwarzschild backgrounds can be found in [8]. There, the Eddington-Finkelstein 
coordinates were shown to be useful for reducing the solution of (1.1) to two elementary 
integrations. In these coordinates, the (background) Schwarzschild metric takes the form 



where v = v(t, r) = t -\-r -\- 2Mln(r/2M — 1) and r = f(r) = r. To investigate the changes 
induced in this geometry by the presence of matter or other fields, wc require a suitable 
ansatz for the complete metric g — g + Sg. For radially symmetric and static source terms 
the most general metric is given by 



< T^i, >G~ G2{G) < Tij,„ >U{1): 



(1.2) 



ds'^ = -(1 - 2M/r)dv^ + 2dvdf + f^dfl'^, 



(2.1) 



ds^ = -e2^(l - 2m/r) dv^ + 26^^ dv dr + dQ'^, 



(2.2) 



3 



where and m are two radially dependent functions to be determined by solving (1.1). This 
reduces to (2.1) in the limit when the source term is shut off. 

The Universe abounds in gauge bosons, and it is therefore of interest to consider in what 
way these particles can alter the spatial geometry surrounding compact massive objects. This 
question can now be investigated by making use of a recent calculation of the renormalized 
electromagnetic stress tensor in a Schwarzschild background. This tensor can be decomposed 
into the two pieces 

< > = < >analyUc +^2(47^^)4 (^.3) 

where < > analytic is given in closed form and AJ^ comes from a numerical evaluation of 
a mode sum [7]. Since this latter contribution is very small in comparison to the analytic 
part, we will not include it in our calculations. The analytic part of the tensor is a simple 
sixth-degree polynomial given by (we drop the angular brackets) [7] 

T/ = — -(-3 -6w- 9w^ - 12w^ + 315w^ - 78w^ + 2A9w^), (2.4) 

* 45(87rM)4^ ^ ' 

Tr = , ,,,, (1 + 2w + 3^2 - 76w^ + 295w^ - 54w^ + 285w^), (2.5) 
45(87rM)4^ ^ ^ ^ 

and 

= T^ = — — -— (1 + 2w + 3w^ + AAw^ - 305w^ + 66w^ - 579^^), (2.6) 
45(87r-M) 

where w = 2M/r. It is straightforward to show this tensor is conserved with respect to 
the background: V^T^ = 0. It is crucial this hold in order that the Bianchi identity 
not be violated. Moreover, one sees by inspection that the renormalized stress tensor is 
asymptotically constant. As r ^ oo, T^^ ^ ciia5f(— 3, 1, 1, 1), which is precisely the 

form of a stress tensor for a thermal ensemble of radiation at a temperature Th — l/(87rM). 
For the Hartle-Hawking thermal state, which represents a blackhole in thermal equilibrium 
with radiation, the polynomial expressions for Tj^ reduce at infinity to that of thermal 
radiation in fiat spacetime. This must be the case in order that an observer at infinity 
measure the correct value of the Hawking temperature, Th- Thus, the asymptotic constancy 
of the stress tensor is a physical requirement. The mass-energy density as seen by a static 
observer is p = — T/, and is negative from r = 2M out to r 5.14M, as shown in Fig. 1. 
As we will see below, this "pocket" of negative energy screens the positive mass M of the 
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singularity from observers located in the neighborhood of the horizon. The existence of this 
negative energy density implies the violation of both the weak and strong energy conditions. 
Moreover, the radial tension {—T^) exceeds the energy density in the range 2M < r ~ 3.2M, 
which results in the breakdown of the dominant energy condition. We mention these points 
as (2.3) provides an explicit example of a physical stress tensor which therefore violates all 
the fundamental energy assumptions of classical general relativity [9]. 

When cast in terms of the coordinates (2.1), these components are given by — T/, 



rpr 

f 



T; and 



rpv 
f 



2M/t)-'[T;-TII (2.7) 

and Tq — are unchanged. The corresponding components of the Einstein tensor which 
derive from (2.2) are 

(2.8) 
(2.9) 



_ -2 dm 

„ _ 2 

r dr 



and 



/-, ^ / \d'4^ dm 
r(l — 2m/r)— — 

dr dr 



(2.10) 



Following [8], we introduce the expansion parameter e = h/M"^. (in units where c = G = 1, 
but h ^ 1). Since the stress tensor (2.3) is known to one-loop, it is 0{h), and consistency 
demands we solve for the metric functions only to the same order. We therefore expand 



e'^ = l + eA(r), m = M[l + eS(r)], 



(2.11) 



which when inserted into (1.1) yield the linearized Einstein equations 

2Me dB 



dr 



2e dA 
r dr 



= 871(1 -2M/r)-i(r;-r/), 



and 



2e 



^^,dA ^dB 
2M)- M' 



87rT;. 



(2.12) 
(2.13) 

(2.14) 



dr dr 

The substitution of (2.12) and (2.13) into (2.14) reproduces the identity (2.7), while the 9 
and components of (1.1) are automatically satisfied since the corresponding constraints 
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VjuT^ = V^T^ are vacuous. Thus, the backreaction problem reduces to solving the two 
equations (2.12) and (2.13). Using (2.4) and (2.5) we integrate (2.12) and (2.13) to find that 
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KB{r) 
and 



-W^ - 13t(;^ + lOSw - ^Iniw) H h ^ H ^ 

3 w 3w'* 



-248 + irS(2M), (2.15) 



KA{r) = ^ 



15 3 1 
3w^ + —w'^ + lOw - 6ln{w) H h 



-32 + X^(2M), (2.16) 



2 ' ' w 2w2. 

where K — 38407r. The gravitational redshift is given by A(r), and the constant of integra- 
tion A{2M) may be fixed after imposing boundary conditions at the radius governing the 
domain of validity of the perturbation (for a detailed discussion of boundary conditions, we 
refer the reader to Sections VII - IX of Reference [8]). However, the remaining discussion 
is independent of this constant, and we do not need to specify it. The function B{r) plays 
the role of a mass-energy function for the gauge fields and is responsible for determining the 
shape of the spatial geometry encoded in the metric (2.2). As discussed in [8], the integra- 
tion constant B{2M) renormalizes the bare mass M via M Mren = M[l + ei?(2M)], and 
can thus be absorbed into the mass parameter. We shall do so with the understanding that 
wherever M appears, it now stands for the "dressed" blackhole mass Mren, ^-nd the constant 
plays no further role in our calculations. The resulting total effective mass-function is 

meff{r) = M[l + eii{r)], (2.17) 

where /i(r) = B{r) — B{2M) is plotted in Fig. 2. This shows how the energy density — T/ 
leads to a negative bosonic mass density in the vicinity of the blackhole horizon. This density 
vanishes at that horizon, achieves its most negative value X/x(r) ^ —120 at r f^i 5.2M 
and passes through zero at r ~ 10. 8M, after which it remains positive definite with the 
asymptotic behavior ~ (r/M)^. To reveal the nature of the modified geometry induced by 
(2.3), we transform (2.2) back to diagonal form with coordinates (t, r, 9, (p) by taking (where 
m stands for me// from here on) 

1^ = 1 and 1^ = e-^(l-2m/r)-\ (2.18) 

which yields the metric 

ds"" = -e2^(l - —)dt^ + (1 - —r'dr^ + r'dn'. (2.19) 
r r 



With the solution (2.19) at hand, we now turn to examine some of its quahtative features. 
A convenient way to visuahze the solution is by "lifting" out the spatial geometry at a fixed 
time, t = constant, and embedding it in a higher dimensional space. The geometry is 
spherically symmetric, so without loss of generality, we confine attention to the equatorial 
slice 9 — 7r/2. The line element of this shce is given by 

ds' = (1 _ dr' + ^2 ^^2^ (2.20) 

r 

and can be embedded in a three-dimensional embedding space. Typically, a single Euclidean 
space would suffice [4,10], but if the slice represents a space of negative curvature, one must 
use a flat Minkowski space (for example, the 3-geometry of Friedmann Universes of negative 
spatial curvature embeds in a 4-dimensional Minkowski space [10]). In either case, the line 
element of the (flat) embedding space is 

dS^ ^±dz^ + dr^ + r^d(f)^, (2.21) 

when written in terms of cylindrical coordinates {z,r,(f)). Identifying the (r, 0) coordinates 
of (2.21) and (2.20) yields the equation for the embedding function 

where z = z{r) decribes a surface of revolution (in the angle (p) about the 2;-axis. For ex- 
ample, in the limit m M (i.e., e — > 0), the left-hand side of (2.22) is positive definite 
(for r > 2M), hence z embeds completely within a Euclidean space. This reproduces the 
standard picture of the Schwarzschild solution with the derivative diverging at the horizon 
r = 2M and going to zero as r — > oo [10]. The square root must be taken in solving for 
z, and both branches meet up at the Schwarzschild horizon, indicating the presence of a 
(non-traversable) wormhole with a throat radius equal to the horizon radius, as expected. 
However, Eq.(2.17) with // ^ C2A* inspection of Fig. 2 indicates that the com- 

bined blackhole- plus-boson mass function m(r) can become negative before reaching the 
Schwarzschild horizon. To embed this region in three dimensions one can use a Minkowski 
embedding space (i.e., choose the minus sign in (2.21)). Outside that region, the solution 
reveals positive definite curvature, and can again be embedded in a Euclidean space. To 
embed all of z at once therefore requires use of a complex embedding space. A somewhat 
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simpler picture containing the same information is had by plotting instead (dz/dr)^, so that 
all regions can be displayed in a single graph. This is done in Fig. 3, which shows a typical 
curve for a solution with a repulsive region. For comparison, we have included the curve 
for the Schwarzschild (e = 0) case. In the former case, the dotted segment represents the 
antigravity domain. There are four length scales associated with this curve. Namely, the 
(renormalized) horizon at r = 2M, the two points where the total mass crosses zero at ri and 
r2, and a cutoff radius rdom which serves to define the domain of validity of the calculation 
(see below). Although the repulsive contribution from fx{r) is too feeble to overcome M for 
a single U(l) boson, the effect becomes pronounced for large numbers of extra abelian gauge 
fields or for large non-abelian groups in the limit when the gauge self-coupling is weak. In 
that hmit, for any gauge group G, we have < T^^, >g~ C2{G) < T^t, >, where C2{G) is 
the second Casimir invariant for G in the adjoint representation, and multiplies the abelian 
tensor given in (2.4-2.6). The corresponding embedding function obtains simply by replacing 
H — > C2{G)iJ, in (2.22). We have calculated the extent of this repulsive domain for a number 
of gauge groups and summarize our results in Table I. To stay within the bounds of per- 
turbation theory, e ~ 1, and this implies a minimum C2 for which the repulsion disappears, 
namely C2 ~ 101. Thus, the smallest unitary, orthogonal or abelian groups which can lead 
to repulsion arc approximately SU(50), SO(50) and [[/(1)]^°°; the corresponding ri and r2 
are given in Table I. We point out that string-inspired models containing large numbers of 
extra abelian gauge fields are not entirely out of the question. Indeed, both Eg x [C/(l)]^^^ 
(C2 = 308) and [C/(l)]^^^ yield anomaly-free superstring theories (in ten-dimensions) and can 
lead to four-dimensional gauge theories that contain more than enough extra gauge fields to 
generate a repulsive domain [11]. Another anomaly-free group relevant for string-inspired 
model building is x Eg, with a C2 = 120 (Table I). For a given C2, there is a minimum e 
which leads to a repulsive region; these are given in Table I. 

To put these calculations in perspective, we should determine the radial domain of valid- 
ity of the metric perturbations leading to the solutions hsted in Table I. That our solution 
must be truncated at some value of the radius can be understood as follows. First, recall 
the stress tensor is asymptotically constant, so it is clear that the corrected metric com- 
puted from (1.1) cannot be asymptotically fiat. However, our calculation yields a solution 
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of the back-reaction (1.1) only to lowest order in h, and must therefore be regarded as a 
strictly perturbative result. The "complete" solution of (1.1) can in principle be obtained 
through an iterative approximation scheme whereby corrections in the metric are fed back 
into the stress tensor. One then recomputes the corrected metric, feeding these higher order 
corrections back into the stress tensor, and so on. At these higher orders, one should take 
into account the quantum fluctuations in the metric itself, by adding an effective graviton 
stress tensor to the right hand side of (1.1). But, to the order we are working here (0(^)), 
it is reasonable to ignore the quantum metric fluctuations when one is computing the semi- 
classical back-reaction caused by the gauge bosons alone to first order. The consequence of 
this approach is that the back-reaction (i.e., the change in the metric) is only meaningful 
in a perturbative sense, that is, the effect of the gauge boson stress tensor is regarded as a 
perturbation of the classical Schwarzschild geometry. In order to regard as a perturba- 
tion there are two length scales which are relevant. One is clearly the mass of the blackhole 
itself. The other scale comes about because, as already pointed out earlier, the corrected 
metric is not asymptotically flat. Thus, in order that the metric corrections be small com- 
pared to the background metric, one cannot consider the perturbations to be of unbounded 
extent. Beyond this cutoff, we can match our solution to an asymptotically flat spacetime. 
Naturally, if (1.1) could be exactly solved, no patching would be expected, and the corrected 
metric would presumably be valid for all r. 

To find this cutoff, we demand that the magnitude of the ratio of the perturbation to 
the background metric be small: 

\^\ = 5<1. (2.23) 
9 

Subtracting (2.1) from (2.2) yields the components of the metric perturbation. 



Sg^^ = 2eC2(G) 



M , , , 2M, / 

—Bir) - 1 U r 

r r 



(2.24) 



Sg,r^eC2{G)A{r). (2.25) 
It is easy to check that both perturbations have the same asymptotic magnitude: 

1^1 ^ !^ (—\^ 
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Inserting this limiting form into (2.23) we find 

The domain radius r^^o^ is defined to be the upper bound for which the metric ratios in 
(2.23) are uniformly small, i.e., \6g\ < \g\ for 2M < r < r^om- To get a fcchng for the 
domain size, consider the case where S — e. Solving for rdom/M from (2.26) shows that 
the domain radius ranges between roughly 8M and 27M (Table I) and we see that all the 
antigravity regions fall comfortably within their respective domains for C2 ~ 1000. This 
gives us confidence that these solutions represent fairly reasonable approximations to the 
back-reaction problem computed to this order. For a given e, as C2{G) increases, the outer 
radius r2 increases but the radius of the domain of validity decreases. Hence, there must 
exist a maximal group for which the perturbative calculation breaks down. This breakdown 
occurs when r2 ~ r^om- This coincidence limit is achieved for groups with C2{G) ~ 1000, 
and the corresponding radius limit is about r2 ~ 8M (Table 1). Of course, beyond the cutoff 
radius, we expect an asymptotically flat spacetime; one can match the perturbed metric to 
a Schwarzschild metric at a radius r2 < f < rdom using the standard junction conditions 
[10]; the mass parameter for the latter geometry is m(f) > (see Fig. 3). That is, 

-9u{r > f) = 1 - (2.27) 
r 

for the tt-componcnt of the metric outside the domain of validity. With this fixed, the 
rr-component will be automatically continuous across the boundary r — f (the angular 
components are vacuously continuous). 
3. Discussion 

By solving for the back-reaction of abelian and non-abelian gauge bosons on a blackhole, 
we have shown that the induced bosonic negative energy density can produce a net negative 
mass distribution (positive hole mass plus negative bosonic mass) extending beyond the event 
horizon. The repulsive force arising from this energy density can dominate the attractive 
pull of the classical blackhole mass for gauge groups with C2(G') ~ 100, and some of the 
string-inspired symmetry groups appear to satisfy this criterion. These results are important 
because they exhibit, already at lowest order, the phenomenon of localized gravitational 
repulsion or "defocussing" , characteristic of traversable wormholes [4]. Of course, we do 
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not expect a traversable wormhole to form at lowest order, but the fact that the requisite 
localized net negative mass-energy is manifest suggests that a wormhole, perhaps dressed by 
field fluctuations, could form (maybe as the endpoint of blackhole evaporation) if the higher 
order back-reactions (1.1) drive the complete solution to pinch-off between the negative 
mass and "normal" regions, but leaving some of the negative mass at the pinch-off site. 
The following comparison may make this point clear. The class of traversable wormholes 
considered in [12] have the following mass distribution 

m{r) ^msheiiS{r -r) + MQ{r -f), (3.1) 

where msheii < and M > 0. Junction conditions are used to match two Schwarzschild 
geometries (of mass M) at the boundary f > 2M. The region containing the blackhole 
(r < f) is excised by hand; thus the wormhole throat is located at f, precisely where the 
total net mass goes negative (note, by the way, the derivative (dz/dr) is finite near the throat 
but discontinuous at r = f, unlike the canonical Schwarzschild case, where it diverges on the 
throat). A glance back at Fig. 2 reveals a qualitative similarity between (3.1) and the mass 
function derived here. As mentioned earlier, we also match on a Schwarzschild geometry at 
for r2 < f < room (of mass m(f) > 0, see Fig. 3). Thus, our mass distribution (2.17) also 
will also include a step function term, m(f) ©(r — f), and the negative mass "well" in Fig. 
2. is qualitatively similar to the the first term in (3.1). 

An analysis of the domain of validity of our calculations shows the location of the re- 
pulsive region lies within this domain for a range of the perturbation parameters. This 
is comforting and suggests the trend revealed by these calculations may persist at higher 
orders. The next iteration step in the back-reaction will be sensitive to the negative energy 
density of the gauge bosons and should be important when the magnitude of the integrated 
mass located between ri and r2, J^^ fi{r)d^r, is of order M. 

The limit e ^ 1 corresponds to the Planck mass limit M Mpi. We see from Table I that 
the radii r2 increase, for a fixed C2, in this limit. The blackhole will of course evaporate, 
but the endpoint of this evaporation could lead to a wormhole (for a large enough C2). 
Although we have not taken into account the effect of quantum metric ffuctuations, it would 
be interesting to include their effect by introducing a renormalized graviton stress tensor on 
the right hand side of (1.1). Such a tensor has in fact been calculated in an arbitrary static 
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vacuum background, but requires "tailoring" before it can be incorporated in the present 
calculation [13]. If we are merely interested in the sign of the graviton's renormalized 
energy density, a shortcut and clue are provided by the fact that primordial gravitons are 
"squeezed" by cosmological expansion [14]. Applying the analysis of [6] implies they will 
have a negative energy density, and so we can hope that the graviton contribution to (1.1) 
enhances the repulsive effect of the gauge bosons. 

Conformally coupled = 1/6) massless scalars also contribute to the dilation effect, but 
the depth of the corresponding pocket of negative energy is roughly one hundred and twenty 
times smaller than the vector boson case and the range is down by a factor of three [15] (The 
back-reaction of conformally coupled scalars on a blackhole has been treated by York [8]). To 
explore the scalar back-reaction case completely, we need to know the scalar stress tensor in a 
Schwarzschild background for arbitrary coupling ^, where is the non-minimal coupling 
of scalars to gravity. To our knowledge, this tensor has not yet been calculated. Scalar 
feedback becomes particularly interesting when looking for wormholes in certain classes of 
higher-derivative gravity theories [16]. In those models, the higher powers of the curvature 
are mapped onto an equivalent scalar sector and the scalar couples minimally (i.e, ^ = 0) to 
the metric. 
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TABLES 



Gauge group: C2{G) e,{emin) n/2M r2/2M rdom/2M 

{6 = e) 

101 (1.) 2.53 2.60 13.4 

120 {Es X Eg) (.85) 2.34 2.81 12.3 

1. 1.73 3.68 

308 {Es X [?7(1)]248) (.34) 2.38 2.75 7.7 

1. 1.15 4.78 

1000 (.11) 1.93 3.38 4.3 

.2 1.30 4.60 

TABLE I. Range of negative mass region and radius of domain of validity vs. gauge group 



FIGURE CAPTIONS 

Fig. 1. Renormalized energy density for an abelian vector boson as a function of ^ 
Fig. 2. The effective bosonic mass-function. 

Fig. 3. Derivative-squared of the embedding function. The dotted hne represents the 
negative mass region between ri and r2 (for e ~ .3 and C2 ~ 300). An asymptotically flat 
geometry (e = 0) can be joined on at a cutoff radius r2 < f < raom- 
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